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^N) ' Abstract. In this brief note, we consider a wave equation that has both 

trapping and a complex potential. For this problem, we prove a uniform bound 
on the energy and a Morawetz (or integrated local energy decay) estimate. 
The equation is a model problem for certain scalar equations appearing in the 

^vj ' Maxwell and linearised Einstein systems on the exterior of a rotating black 

(N ■ 



^ 



hole. 



■*^ ■ 1. Introduction 

r^^ ' We consider the Cauchy problem 

^ • {-d^ + dl + V^(A^ -N)+ leW) u^O, (1) 

u{0,x,uj) ^tl)Q{x,uj), dtu{0,x,uj) = tlJi{x,Lo), (2) 

on (t, x,uj) G A/ = K X M X S*^ with smooth, compactly supported initial data. Here 
^ ' u is a complex function u = v + iw, 

i>: 1 

\n . X + i 

^^ ' VF is a smooth, real-valued, compactly supported function which is nonvanishing at 

r«^^ , X = and uniformly bounded by 1, and e > is a small parameter. Finally, A^ is 

(^ ' the Laplacian in the angular variables and A'^ is a number chosen to be sufficiently 

large to allow us to avoid certain technical problems. 

The equation ([1]) has both trapping, which occurs at x = 0, and a complex 
potential, which does not vanish at the trapped set. The interaction of these creates 
problems, which appear to frustrate the use of energy and Morawetz estimates at 
the classical level. By adapting known, pseudodifFerential methods, we show how 



C^ ' to overcome these problems. We now state our main result in terms of the energy 



m - 1 f 



\dtuf + |9,u|' + V{\Vujuf + N\uf) dxd'uj, 

<S2 



Theorem 1. There is a constant C such that, if ipo and ipi are such that E{0) is 
finite then 

Vi e M : E{t) < CE{0), (3a) 

\d.u\- , X- ( iV.uP , \dM' ^ , \u? ^^^^^.^ < ^^^0)^ ^3^) 

M 



x'^ + l i-^x^ll + kP x'^ + l) l-t-lxl' 



/ M^ dtdxd^c. < Ci?(0). (3c) 

Jm 1 + \xV 



Date: July 25, 2011 File: ModelProhlem.tex. 

1 



2 L. ANDERSSON, P. BLUE, AND J. -P. NICOLAS 

Since the equation ([T]) has t independent coefficients, one might naively think 
that Noether's theorem provides a positive conserved energy. However, for the 
Lagrangian Ci[u,du] = —{dtuY + (d^u)^ + Viy^jU ■ V^u + Nu^) — ieWu'^, which 
has the wave equation ([1]) as its Euler-Lagrange equation, the conserved quantity 
associated to the time translation symmetry is indefinite, being approximately the 
energy of the real component of u minus the energy of the imaginary component 
(plus e times a term involving Wvw). On the other hand, a Lagrangian of the form 
C2[u,du] = — |9tup + |9a;Wp + V^(|VijUp + iV|up), which corresponds to the energy 
expression considered above, does not yield equation ([T]) as its Euler-Lagrange 
equation. 

The wave equation ([T]) is a model for equations arising in the study of the Maxwell 
and linearised Einstein equations outside a Kerr black hole. The Kerr black holes 
are a family of Lorentzian manifolds arising in general relativity, and they are 
characterised by a mass parameter M and an angular momentum parameter a. 
Black holes are believed to be the enormously massive objects at the center of 
most galaxies. The case \a\ < M is the physically relevant. The case a = is the 
Schwarzschild class of black holes. 

It is expected that every uncharged black hole will asymptotically approach a 
Kerr solution under the dynamics generated by the Einstein equations of general 
relativity. The wave. Maxwell, and linearised Einstein equations on a fixed Kerr 
geometry are a sequence of increasingly accurate models for these dynamics. By 
projecting on a null tetrad, the Maxwell and linearised Einstein fields can be decom- 
posed into sets of complex scalars, the Newman-Penrose (NP) scalars [22| , [23 l IT4] . 
It is well-known that the NP scalars with extreme spin weights satisfy decoupled 
wave equations, known as the Teukolsky equations, and that the solutions to these 
reduced equations can be used to reconstruct the full system [29] . 

For the Maxwell field on the Kerr background, the spin weight NP scalar can be 
treated in the same way, and the resulting equation is known as the Fackerell-Ipser 
equation [llj . For linearized gravity on the Schwarzschild background, it is also 
well known that the imaginary part of the spin weight NP scalar is governed by a 
wave equation, the Regge- Wheeler equation [26^, 2^ ■ The corresponding equation 
for the real part is more complicated, cf. [33l[2T], see also [I]. 

It was recently shown p] that in the general (|a| < M) Kerr case, by imposing a 
gauge condition related to the wave coordinates gauge, the equation for both the real 
and imaginary parts of the spin weight NP scalar of the linearized gravitational 
field may be put in a form analogous to the Regge- Wheeler and Fackerell-Ipser 
equations. Explicitely (in the Kerr spacetime with signature — \--\ — h, working in 
Boyer-Lindquist coordinates) these take the form 



M 
^"^" + 2s2— — — -^)u = 0, (4) 



where s = corresponds to the free scalar wave equation, s = 1 corresponds to the 
Maxwell (Fackerell-Ipser) case, while s — 2 corresponds to the linearized gravity 
(generalized Regge- Wheeler) case. In particular, for the a ^ cases, the analogues 
of the Regge- Wheeler equations have complex potentials, with the imaginary part 
depending continuously on a. 
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For the wave equation in the Schwarzschild case, the use of the energy estimate 
[31], hke (j3a| with C = 1, and Morawetz estimatcfUl are weh estabhshed [171 El 
m ini IB]- In the Morawetz estimate (j3bl) . there is a loss of control of time and 
angular derivatives near a; = 0, in the sense that the integrand cannot control 
\x\P{\u\i'\dtu\^-'i* + |u|9"|V^u|2-«-) with both p = and either g^ = or g„ = 0. 
The presence of trapping makes some loss unavoidable [25] . By applying "angular 
modulation" and "phase space analysis" , the range for the angular parameter q^^ 
can be refined to p = and g > [5]. This type of refinement is crucial in the 
current paper. Alternatively, certain pseudodifferential operators have been used 
to obtain refinements near x = 0, to p > and qt — Qi^ = [201 . 

For the wave equation in the general (|a| < M) Kerr case, it is possible to apply 
Fourier transforms first in the (j) and t variable^ and then the remaining variables. 
The individual (f> modes decay pointwise |il2;. Although the problem has a time- 
translation symmetry, because the generator of time translations fails to be a time- 
like vector with respect to the Lorentzian inner product of the Kerr geometry, there 
is no positive, conserved energy. A major advance was the proof that, in the slowly 
rotating case \a\ <^ M, there is a uniform energy bound, like estimate ([3a|) . The 
first proof used an estimate similar to (|3b| , but with additional restrictions on the 
support of the Fourier transform [7 . Independent work [55] established estimates 
similar to ([5a|) and pbp . but with no restriction on the Fourier support, and there 
were subsequent pseudodifferential refinements |30]. Also, the first two authors 
have proved similar results using methods which require two additional levels of 
regularity but which completely avoid the use of Fourier transforms. Morawetz 
estimates and refinements are a crucial step in proving pointwise decay estimates 
[SI [5] m m [m [H] and Strlchartz estimates [IHllSSl, including the long-conjectured, 
inverse-cubic. Price law [l6l[27] . 

The study of the Maxwell and linearised Einstein systems in the Kerr geometry 
is still in its infancy. For the general Kerr case, a certain transformed, separated 
version of the Teukolsky system has no exponentially growing modes [32j . In the 
Schwarzschild case, the (p modes of the Teukolsky equation decay pointwise [13| . 
Recently, improved decay estimates for the Regge- Wheeler type equation ([4]) on the 
Schwarzschild background, giving decay rates of t^'^, t^** and t^^, respectively, for 
s = 0, 1, 2, have been proved [10]. The Regge- Wheeler equation has been used with 
the full system to prove energy and Morawetz (and pointwise decay) estimates for 
the Maxwell system [5] and, more recently, under assumptions on the asymptotic 
behaviour, for the full (not merely linearised) Einstein equation |I5\ 

For the spin-weight equations arising from the Maxwell and linearised Einstein 
equations with a ^ 0, the presence of complex potentials in the reduced equations 
l[4| prevents the existence of a positive conserved energy and means that an unre- 
fined Morawetz estimate, such as (j3b|) . is insufficient to control the growth of the 
energy, which is why we also prove estimate (I3cl) . In the Kerr geometry, there is 
no positive, conserved energy because the generator of the time-translation sym- 
metry fails to be timelike everywhere. In contrast, for the equation ([T]), there is no 
positive, conserved energy because the complex potential prevents any variational 
approach from providing an energy-momentum tensor that satisfies the dominant 



These are also called integrated local energy decay estimates. 

Here is the azimuthal angle, which would be one component of u) in the notation of this 
paper. 
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energy condition. The method of this paper combines a Fourier-transform-in-time 
technique (as in [71 [28]) with a "modulation" (or Fourier- rescaUng) technique (from 

As is common, C will be used to denote a constant which may vary from line to 
line, but which is independent of the choice of u or T. The notation ^ < _B is used 
to denote that there is some C such that A < CB, with C independent of u and 
T, and similarly for >. 



2. A PRELIMINARY ENERGY ESTIMATE 

We derive an estimate for an energy for the wave equation ([1} by integrating 
by parts against dtu and following the standard procedure for getting an energy 
estimate: 

= Re {{dtu){-df + d^ + V{A^ - N) + ieW)u) 

+ dxRc{{dtu)dxu) - -dt\dxu\'^ 

+ V^ • Re{{dtu)V^u) - ^dt {VilV^uf + N\uf)) 
— eWlm{{dtu)u) . 

Introducing an energy which we denote by 

Eit) = i / \dtuf + |9,up + V{\V^u\^ + N\u\^) da;d2w, 

^ J{t}xRxS^ 

assuming that u decays sufficiently rapidly as |a:| — ?> oo, and integrating the previous 
formula over a region [f i , ^2] x K x S*^ , we find 

E{t2)-E{ti)= I -eWlm{{dtu)u)dtdxd^uj. (5) 

"'[ti,t2]xRxS2 

In particular, note that the energy fails to be conserved and that an estimate of the 
form (j3bp would be insufficient to control the right-hand side. There is, however, a 
trivial exponential bound: 

<t2-tl). 



E{t2) < e'^*'-''>E{ti) 



3. The Morawetz estimate 



Following the standard procedure for investigating the wave equation, we derive 
a Morawetz estimate by multiplying the wave equation by {f {x)dxU + q{x)u) , where 
/ and q are real- valued functions. 

In performing this calculation, it is useful to observe that 



q' (x)Re{udxu) = dx i —uu j - -q'uu. 
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Using this and applying the product rule term-by-term, one finds 

Re( (fd^u + qu) {^dfu + d'^u + V{A^ - N)u + ieWu) ) 
= dtpt + dxPx + yujPu 

-Ir + q) \dM' - (Ir + q) \dxv 



1 ./ ^ T. , 1 ,r^ TM \ ,W ,2 



^f -1 y+2^id^v)]\^^^\ 



'■" ^ \u? 



efWlm{{dxu)u), (6) 



where 



Pt =ptif,q;u) = -Rc{{ f {dxu) + qu){dtu)), 

Px = pM,q, «) - lmu\' + lf\dxu\' - ^fV\\/^u\' 

+ qReiudxu)-^iNfV + q')\u\'', 
Pc - pM, g; u) - fVRe{{dxu){V^u)) + qVRe{uV^u). 

We take / = -arctan(a;), for which /' = -{x'^ + 1)-^ = -V, f" = 2x{x'^ + 1)-'^, 
and /'" = -2(3a;2 - l){x^ + 1)-^. We take q = /'/2 + (5(1 + x^)-^ arctan(a;)2 for 
some sufiieiently small S. 

We use the notation 

Efd^+qit) = / Re{f{dxu)dtu) + Re{qu{dtu)) da;d2w, 

J{t}xRxS^ 

and observe that, by a simple Cauchy-Schwarz argument, there is the estimate 
\Efo^+,\ < CE. 

Observing that the left-hand side of ^ vanishes, we have 



= dtPt + dxPx + ^u:Pu. 

-^ "^^"""^7 \dtu\^ + — ^(l-(5arctan(a;)2)|a^up 



,arctan(x)^ , „ ,9 1 

1 + X'' 1 + X'' 



l^uU 



2 



X arctan(.T) — S arctan(x) 

/a; arctanfa;) — (5arctan(a;)^ \ 1 „\ , ,, 
N \' ,,, ^-^ + 7:1 H 



(l + a;2)2 ; ' 2 

— efWIm{{dxu)u) 

Taking e sufficiently small, N sufficiently large, and 6 sufficiently small, the factors 
in front of |i9a;Mp and jup are nonnegative and one can dominate the term involving 
W using these two terms. (These estimates are uniform, in the sense that, if the 
estimate holds for choices of eo , Nq, and Sq , then it remains valid for e < eq, N = Nq, 
and S = Sq.) Thus, by integrating over a time-space slab M[t^ j^] = [^ij ^2] x M x S^, 
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one can conclude that there is a constant C such that 

E{t2) + E{h) 



^ / V-T + I arctan(x)|" f-^ViT + V-^ + , ,, dMxd^w. (7) 

4. PSEUDODIFFERENTIAL REFINEMENTS 

4.1. The wave equation for an approximate solution. We define a smooth 
characteristic function of an interval [a, h] to be a function which is identically 1 on 
[a, h\, which is supported on [a — 1, 6 + 1], and which is monotonic on each of the 
intervals [a — 1, a] and [&, 6 + 1]. A smooth characteristic function of a collection of 
intervals, each of which are separated by distance at least two, is defined to be the 
sum of the smooth characteristic functions of each interval. 

Let T > be a large constant. (Here large means larger than — log|e| and 
2.) Let xi be a smooth characteristic function on [0,T], and let X2 be a smooth 
characteristic function of [—1, 0] U [T, T + 1]. Let x\x\<2 be a smooth characteristic 
function of [—1,1]. We will use xij X2: and x\x\<2 to denote Xi(^)j X2{i), and 
X\x\<2{x) respectively. 

Since xi is smooth, there is a uniform bound on its derivative and second deriv- 
ative, each of which are supported on [0, 1] U [T, T + 1], so that there is a constant 
C such that l^tXil + \d'ixi\ < C'X2- 

The functions 

Wl = XlX\x\<2'U; 
U2 = X2X\x\<2U, 
U3 = XlU, 

satisfies the equation 

{-dl +dl + T/(A^ -N)+ leW) ui ^ F{u2, Vu2,t, x) + G{u3,\7u3, t, x), (8) 

where 

F{u2,\7u2,t,x) = -2{dtXi){dtU2) - {dtXi)u2 

G{u3, Vm3, t, x) = 2{dxX\x\<2){dxU3) + {d'^X\x\<2)u3- 

Since all functions oft in this equation are smooth and supported in t G [— 2,T + 2], 
they are Schwartz class in i, so we may take the Fourier transform in t and remain 
in the Schwartz class. We will use ^ to denote the Fourier transform in i, and r 
for the argument of such functions. We will typically use the word "functions" to 
describe u, ui, U2, and U3 and the words "Fourier transforms" to describe their 
Fourier transforms. We will use L^ to denote i^(da;dxd/:) for functions and to 
denote i^(da;dxdT) for Fourier transforms. We will use || • || for || • ||^2 unless 
otherwise specified. 

We introduce the following space-time integrals 

I{T)^ I I I x^\dtu\^ + \d,^u\'^ + \u\'^ALudxAt, 
J -2 J-2Js^ 



J{T) = ( Irl^/^l^pdcjdxdT. 

jRxRxS^ 
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The dependence of J upon T is through the smooth cut-off xi in wi. Typicahy, the 
argument T wiU be clear from context and wiU be omitted. From the Morawetz 
estimate ([7]) and the exponential bound on the energy, it follows that / < E{T) + 

We now aim to prove a Morawetz estimate using the Fourier transform. We take 

/ ~ — arctan(|r|"a;), 
_ /' _ 1 |t|" 



2 21 



|2a^2 ■ 



with a G [0, 1/2]. We multiply the Fourier transform of equation ([8]) by {fdx + q)ui, 
and integrate the real part over M x R x 5^. This integral is convergent because 
all the functions are compactly supported in time, so the Fourier transforms are 
Schwartz class. 



4.2. Controlling the terms arising from the cut-off. We consider first the 
integral arising from the right-hand side of ([5]). This is 

Re({{fdx + q)^)(F + G))dujdxdT<\\{fdx + q)^i\\\\F + G\\. 

The terms on the right can be estimated by 

||(/5, + g)^||< 11/5,^11 + ||g^||, 
11/5.^11 < 115.^11 </'/', 

ik^ii < iikr^ii < iRii + iiirr/^^ii < /1/2 + Jl/^ 

and 

||i? + G||<||F|| + ||G||. 
Because G is supported only for t e [—1, T + 1] and x e [—2, 2], we have 

\\G\\<I'/'. 

Similarly, because F is supported only for t G [—1,0] U [T, r+ 1] and x £ [—2,2], 
we have that at each instant in t, the function F is bounded in L'^{dxduj) by either 
GE{Oy/^ or GE{TY/'^ . Since we are considering two intervals in t of length 1, we 
have 

ll^ll - ll^ll ^ c(£;(o)i/2 + e{tY'^). 

Thus, the terms on the right-hand side of the Fourier transform of ([8]) are bounded 

by 

/ Re {{fdx + q)^) (f + g) dujdxdr 

< c(£;(o)i/2 + £:(T)i/2 _^ ji/2)(£;(o)i/2 _^ Eirf^). (9) 
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4.3. The Morawetz estimate for the approximate solution. If we multiply 
the left-hand side of the Fourier transform of the wave equation ([8]) by {fdx + q)ui 
and take the real part, then we have the analogue of ^ 

Re( (fdxui + qSi) {t^ui + d^ui + ^(A,^ - N)1Ii + ieWu{) ) 

= dxPx + y^puj 

+ ( Q/'- -7)^^ + ^/(5.1^)) I v.^P 

A^ ( (^/' - 9) ^ + Ifid.v)^ + ^g") RP 

-efWIm{{dx€i)^), (10) 

where 

Px =P:,(/,g,^) = -/|9t^P + -/|9x^P - -/V^lV^M^p 

+ qRet^A^) ~ ^{NfV + q')\^\\ 
Pu =Puj{f,q;ui) = fVRe{{dxul){V^ii^)) + qVRe{ul\/^ul). 

Note that there is no pt term because, for Fourier transforms, the analogue of the 
product rule is simply — iruiMi — uiirui. 

When this equality is integrated over a space-time slab, the px and p^j terms 
integrate to zero, and the remaining terms are all non-negative except for those 
arising from q" and from W. The integral of the term involving W is bounded by 
/. 

We now consider the term involving q": 

2^'^' '' (l + |r|2":E2)3l"il • 

From the positivity of the remaining terms and the bound ([9]) of the terms coming 
from the right-hand side of the wave equation ([8]) for ui, we have that 

/ |rp° ,^ ~ f^rt, |^|^dc.da:dr 

JrxRxS- (l + |r|2"x2)3l 

< c(£;(o)i/2 ^ £:(r)i/2 _^ ji/2)(£;(o)i/2 ^ E{Ty/^). 

This can be combined with an additional factor of MI, where M is a large constant. 
(700 is sufhcicnt.) The integral / dominates the integral of (|t|2 -|- l)a;2|u^|2 and is 
bounded by C{E{T) + E{Q)). Thus, we have 



/ 



<C{E{Oy/^ + E{Ty/^ + ji/2)(^(0)i/2 ^ EiTf/^). 
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By considering the two cases |t|"|x| < M ^1"^ and |r| 
that if 2 — 2a = 3q (i.e. a = 2/5), then 



> M ^/^, one can see 



.|3a 



1 



\2a^2\3 



x'^y 



(1 + |T 

and therefore we find 

C{E{Oy/'^ + EiTf/^ + ji/2)(^(0)i/2 + E{TY/^) 
> [ /|_|3. l-3|rp"a:^ 

which imphes 



+ M(|t|2 + l)x^ X|.|<2 > C|t|3"x|.|<2 



Af(h 



J <C{E{T)+E{Q)). 



l)x^ l^l-'dwdxdT > J, 



(11) 



4.4. Closing the energy estimate. It is now possible to estimate the integral 
on the right-hand side of the energy estimate ([5]). For |a;| > 1, the right-hand 
side (using the compact support of W and the Cauchy-Schwarz estimate) can be 
dominated by / < C{E{T) A- E{Q)). For |x| < 1, we would like to dominate the 
integral over M x R x S*^ of \WlTn.(uidtUi)\ by the integral J. However, this is not 
entirely correct, because in J there is a contribution arising from the support of u 
in the region i G [— 1, 0] U [T, T -|- 1]. The error in this approximation is bounded 
by C{E{T) + E{0)). Thus, we have 



E{T) - E{0) < Ce E{T) + E{0) 



Wlin{uidtUi)dujdxdt 



We can also take the Fourier transform to obtain an estimate by 



Ce(E{T) + E{0) 



xRxS2 



W lT[i(iiidtUi)dujdxdt 



< 



Ce (e{T) + E{0) + f \t\ |^pdwd.Tdi 

V JRxExS^ 



The integrand is now controlled by / + J, which, from estimates ([7]) and (fTTj) , we 
know can be estimated also by the sum of the initial and final energies. This leaves 
the estimate 

E{T) ~ E{0) < Ce{E{T) + E{Q)). 

By taking e sufficiently small relative to the constant, we obtain a uniform bound 
on the energy 

E{T) < CE{0). 

We note that, since all the constants were independent of T, the estimate holds 
uniformly in T. This proves the first statement (j3a[) . in theorem [TJ Combining this 
with estimate ([7]) (and estimating x^{l 



< 



arctan(a;) ) gives the second, (|3b|) . 
Finally, the arguments of this section and the bound on / + J, from estimates ([7]) 
and (fTTj). give the third result (pc)) . 



Remark 2. Using the method given in [5], the stronger Morawetz estimate (|11|) 
can be improved to control the integral of |rp~^|up for any £ > 0. Because of the 
presence of trapping, it is not possible to improve this to \t\'^\u\'^. 
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